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FIG. 1. Direct laser writing in nanodiamond photoresist, a, Sketch of the direct laser 
writing process. A femtosecond laser beam is focussed into the photoresist in order to polymerize 
well defined 3D structures, b, Scanning electron micrograph of such a structure after development 
containing several key photonic elements, such as waveguides, couplers and microdisc resonators. 
Scale bar is 5 \xm. 

A fully integrated quantum optical technology requires active quantum sys- 
tems incorporated into resonant optical microstructures and inter-connected in 
three dimensions via photonic wires. Nitrogen vacancy- cent res (NV-centres) 
in diamond which are excellent photostable room temperature single-photon 
emitters 1 are ideal candidates for that purposePS. Extensive research efforts 
to couple NV-centres to photonic structures such as optical microresonatord^El, 
microcavitied^, and waveguide^^ have been pursued. Strategies for integra- 
tion range from top-down fabrication via etching of diamond membraned^to so- 
phisticated bottom- up assembly of hybrid structures using diamond nanocrystaLpSED 
where the latter approach allows for deterministic coupling. Recently, another 
approach based on the incorporation of nanodiamonds in soft glass optical fibres 
via a melting process has been introduced 12 . Here, we utilize two-photon direct 
laser writing (DLW) to fabricate fully three-dimensional (3D) structures from a 
photoresist mixed with a solution of nanodiamonds containing NV-centres. For 
the first time, this approach facilitates building integrated 3D quantum photonic 
elements of nearly arbitrary shapes. 

On-chip photonic circuits working at the single-quantum level play an important role 
for future quantum information processing^. Several approaches to assemble such quan- 
tum photonic devices from different fundamental photonic entities have been pursued^. 
For example, by using self- assembled quantum dots in semiconductor membranes, one can 
exploit the full power of semiconductor nanofabrication technology, and sophisticated struc- 
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tures have been demonstrated 13 , even with on-demand coupling architecture^. However, 
this approach is limited to two-dimensional structures, except for few results obtained by 
extremely challenging manual membrane-stacking 16 . 

Another easy and low-cost way of fabricating photonic structures is optical lithography 
via direct laser writing (DLW)^^ where a tightly focussed femtosecond laser beam is used 
to expose a photoresist (see Figure [jji). The use of multi-photon absorption enables a 
sequential 3D exposure by scanning the sample or the focus of the laser. For common 
negative-tone photoresists, unexposed parts are removed during a development step and the 
3D polymer structures remain (see Figure [TJd). DLW is well known for the fabrication of 
photonic crystals 18 or other photonic elements like resonato nj 19 | 2Q | and waveguide^. In order 
to functionalise the structures with optically active material, fluorescent dyed, quantum 
dotd and metal nanoparticled have been incorporated. However, until today there has 
been no 3D structure operating at the fundamental quantum level with single photons from 
single emitters being collected and routed. Moreover, no combinations of multiple optical 
elements (different resonators, couplers, waveguides) have been demonstrated. The reason 
is the lack of photostable quantum emitters which are compatible with the DLW process 
while still preserving the possibility for high-quality DLW fabrication. 

Here, we overcome these obstacles by mixing small amounts of nanodiamonds from solu- 
tion directly into a known photoresist!^ (See Methods). This special photoresist allows for 
the fabrication of transparent optical elements made from an acrylate polymer containing 
nanodiamonds. The nanodiamonds contain NV-centres which are photostable even after 
exposure by the DLW laser. After the DLW process they serve as integrated single-photon 
sources. It is especially noteworthy that our technique is not limited to NV-centres in dia- 
mond nanocrystals. It can be applied to any photostable single-photon emitter in nanocrys- 
talline material, e.g. other diamond defect centres, such as Silicon vacancy-centred, or 
other deep defect centres in large band gap material. 

In the following, we demonstrate the fabrication, functionality, and interaction of key 
building blocks for quantum photonic circuitry. These are interconnects (i.e., waveguides) 
and functional elements like resonators and emitters. We measure single-photon emission 
from a diamond nanocrystal in a waveguide-coupled resonator which demonstrates the func- 
tional combination of all these building blocks and paves the way to more complex devices. 

As a first component for integrated photonic circuits, disc resonators with a disc diameter 
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FIG. 2. Mode measurements, a, A tunable laser is coupled to a disc resonator via a tapered 
fibre. Upon sweeping the laser frequency, the transmitted light is modulated by the modes of the 
resonator, b, Resonator mode with a quality factor Q of 1.2 x 10 4 . c, Scan over many modes with 
a free spectral range of 6.5 nm indicated by the red arrows. 

of 20 |xm and a disc thickness of approximately 1.2 |xm on a stem with diameter 10 |Lim 
were fabricated using the photoresist functionalised with nanocrystals (see Methods). The 
resonator supports whispering gallery modes which are confined in the disc close to the 
rim with high quality factors (Q-factors). We measured the Q-factor using a tunable laser 
coupled into the resonator via a tapered fibr^S We found Q-factors as large as 10 4 at a 
wavelength of around 770 nm (see Figure [2] and Methods). The observed free spectral range 
of 6.5 nm of the resonator corresponds well to the expected value of 6.4 nm for a disc with 
diameter 20 \xm derived from geometrical considerations assuming an index of refraction of 
1.5. 

Next, a homebuilt confocal microscope (see Methods) was used to raster scan the fab- 
ricated resonator discs and to identify single NV-centres by measuring the second-order 
autocorrelation function g^(r) in a Hanbury Brown and Twiss (HBT) interferometer. In 
parallel, a grating spectrometer was used to resolve the emission spectra of the individual 
emitters. Figure [3)3 shows a confocal scan of a resonator where fluorescent defects can be 
identified as bright spots on the resonator. Encircled is a spot on the resonator's outer rim 
where coupling to the disc's whispering gallery modes is expected. Figure [3]l displays an 
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autocorrelation measurement of fluorescence collected from that spot. A clear antibunching 
dip is observed reaching a value of fl^(O) = 0.31 as deduced from the fit shown as red curve. 
No background correction was applied to these data. This measurement shows that the 
bright spot indeed corresponds to a single NV-centre in the disc resonator. 

The fluorescence spectrum of an NV-centre at room temperature is typically broadened 
over 200 nm from approx. 600 nm to 800 nm by phonon sidebands. This corresponds well 
to the measured emission spectrum shown in Figure [3J3. The peak at approximately 630 nm 
stems from fluorescence of the photoresist and bleaches after long excitation. For confocal 
scans and correlation measurements this background light was suppressed by spectral filter- 
ing. In addition, there is a fine modulation of the spectrum due to the cavity resonances. It 
is attributed to photons that were initially emitted into the resonator modes and afterwards 
scattered out by the diamond nanocrystal. Hence, these photons are detected in addition 
to the flat unmodulated spectrum emitted directly out of the disc. 

As a second photonic element, arc waveguides from the same material with a width of 
1.8 |im and a length of 40 \im) were built. The waveguides were fabricated close to the disc 
resonator in order to enable evanescent coupling to the resonator's modes (see Figure |4^l). 
Their ability to guide light within a broad spectral range is demonstrated by illuminating one 
end of the waveguide with light from a halogen lamp as seen in Figures^ and c. The blue 
and red circular spots correspond to the input and output of the arc waveguide, respectively. 

Coupling via the waveguides to the resonator modes was confirmed by shining an ex- 
citation laser (wavelength of 532 nm) into one port of the waveguide and, after blocking 
excitation light by a dichroic mirror, analysing the light coming out of the other port. Ob- 
viously, the 532 nm-laser excites background fluorescence in the resonator, which is then 
coupled back into the waveguide. The measured fluorescence spectrum in Figure [4^ thus 
shows the characteristic modulation by resonator modes and proves a good coupling between 
guided modes in the arc waveguide and confined modes in the resonator (see Supplementary 
Information) . 

By moving the excitation spot of a second objective over a waveguide, individual nan- 
odiamonds inside can be addressed and identified as spots with higher fluorescence (see 
Figure [i]). The intensity cross correlation between the waveguide's two outputs shows a 
clear antibunching behaviour. This result highlights a way of coupling single emitters to 
thin waveguides, which can possess very high collection efficiencies, as demonstrated re- 
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FIG. 3. DLW resonator containing single NV-centres. A diamond nanocrystal containing a 
single NV-centre is coupled to whispering gallery modes of a DLW disc resonator, a, Measurement 
scheme. Detection and excitation take place at the same point in a confocal configuration, b, 
Scanning confocal image of the resonator disc. The circle indicates a bright spot identified as 
single NV-centre. Its fluorescence is analysed in c and d. Scale bar is 5 \im. c, Spectrum of the 
collected fluorescence. The resonator modes are seen on the broad NV-centre phonon sidebands. 
The peak at 630 nm stems from the photoresist and can be bleached over time, d, Autocorrelation 
function of the fluorescence from the NV-centre. A clear antibunching behaviour can be seen. The 
red curve is a fit to the data. 

cently using tapered optical fibred^S. Here, we show a direct on-chip integration of such 
devices with both waveguide ends accessible. Detection at the waveguide's ends can yield 
count rates of 29, 000 /s after passing a 50/50 beam splitter and spatial and spectral filter- 
ing. The autocorrelation function at zero time delay of the fluorescence collected in this 
way is g^(0) = 0.37 showing the single-photon character of the light (see Figure When 
accounting for the 50/50 beam splitter, at high excitation powers the number of background 
corrected counts from the defect centre is measured to be 130, 000 /s. 

Finally, we verified active quantum functionality of a coupled 3D system, i.e., single 
photon emission to resonator modes, as well as subsequent collection and routing via a 
coupled waveguide. In order to do so, raster scans of the excitation spot over the resonator 
disc were performed while the photons from one of the waveguide's ends are collected (see 
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FIG. 4. Arc waveguide coupled to resonator, a, Scanning electron micrograph of a 3D 
arc waveguide /resonator assembly. Length of the scale bar is 5 \xm. b, The waveguide resonator 
assembly as seen through the cover slip. The upper half of the image is illuminated from below 
by white light. Efficient light guiding through the arc waveguide can be seen. Scale bar is 10 \im. 
c, Sketch of the region shown in b with the illuminated area marked yellow, d, Sketch how the 
spectrum in e was acquired, e, Spectrum of the fluorescence detected at one of the waveguide's 
ports while exciting through the other. Resonator modes can be clearly identified, f, Experimental 
configuration for the autocorrelation measurement in g. A second objective is used for excitation, 
g, Measured autocorrelation function of a NV-centre inside a waveguide. It is acquired at both 
output ports of the waveguide in a cross-correlation configuration. 

Figure [5^i). When the excitation laser hits the outer rim of the resonator, the created 
fluorescence light is coupled to the resonator modes and subsequently out to the waveguide, 
as can bee seen in Figure [5)3. Therefore, the outer rim of the resonator appears bright in the 
raster image. Bright spots indicate possible positions of NV-centres. To identify NV-centres 
via their non-classical light emission, an HBT setup, with the laser-written arc waveguide 
serving as beam splitter, is used. Figure shows this correlation measurement for the 
NV-centre encircled in Figure [5^t. The data shown were acquired with a pulsed laser at 
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FIG. 5. Waveguide coupling of a single NV-centre inside a resonator, a, Sketch of the 
experimental configuration. The excitation spot is scanned over the resonator disc. Photons are 
detected at both waveguide outputs simultaneously, b, Photon counts collected at one end of the 
waveguide while scanning the excitation spot with a second objective. The circle highlights the 
position of a single NV-centre. Shape distortions are due to non closed-loop piezo-scanning. Scale 
bar is 5 \xm. c, Cross correlation (background corrected) between the waveguide output ports when 
the NV-centre marked by a circle in b is excited. The red curve is a fit to the data. 



40 MHz repetition rate and background corrected (see Supplementary Information). The 
value of g( 2 \0) = 0.27 indicates that the main contribution to the signal stems from a single 
NV-centre. This clearly indicates that the NV-centre emits single photons into resonator 
modes, which are then coupled out by the arc waveguide. With a single source in a resonator 
coupled to a waveguide, this device represents a key integrated 3D quantum photonic circuit. 

In conclusion, we have shown an approach which allows for the direct incorporation of 
single quantum emitters into true 3D photonic structures of nearly arbitrary shape. The one 
step fabrication scheme is easy to implement and requires no elaborate clean-room environ- 
ment. We have pointed out that our method is not limited to NV-centres in nanodiamond, 
but can be extended to any other stable quantum emitter embedded in nanocrystals. The 
flexibility of DLW enables immediate scaling-up to more complex structures (see Figure [TJd) 
where single photons are collected, sent through beam splitters, interferometers, or other 
optical elements inter-connected by waveguides in a 3D architecture. Tuning of individual 
components, e.g. the resonator, is possible as well (see Supplementary Information). With 
(partial) metallisation of the structures embedding single photon emitters in 3D plasmonic 
structures or in metamaterials 30 can be envisioned. A next challenging, but rewarding step 
is to perform DLW of perfectly aligned 3D structures around emitters of known position. 
This would enable on-demand fabrication of arbitrary 3D quantum photonic architectures. 
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Methods 



Direct laser writing. The DLW setup consists of a Ti:sapphire oscillator (Spectra- 
Physics MaiTai HP) delivering 100 fs pulses centred around 810 nm. The beam is focussed 
through an oil-immersion lens (Leica HCX PL APO 100x/1.40-0.70 OIL CS). The sample 
is moved with a 3D piezo stage (PI P-527.3CL). The writing was performed with 6 mW av- 
erage writing power and at 50 [im/s scan velocity. The samples were developed in isopropyl 
alcohol for 30 minutes, afterwards in acetone for 30 minutes, and finally dried super-critically 
(Leica CPD030). 

The new photoresist is based on the monomer pentaerythritol tetraacrylate (PETTA) 
which contains 350 ppm monomethyl ether hydroquinone as inhibitor. Next, 0.25 % wt 
of the photoinitiator 7-diethylamino-3-thenoylcoumarin were added. Last, 2 % wt of an 
ethanol-based nanodiamond suspension were added and stirred overnight. Nanodiamonds 
were of type lb and had a median diameter of 25 nm (Microdiamant AG). Spatial dimensions 
of the fabricated structures were extracted from scanning electron micrographs. 

Confocal microscope. The homebuilt confocal microscope is based on the body of a 
Zeiss Axiovert 200. Either a frequency doubled Nd:YAG continuous wave laser or a pulsed 
laser (LDH-P-FA-530, PicoQuant) were used to create the excitation light around 532 nm. A 
large-numerical-aperture air microscope objective mounted on a ID piezo stage (MIPOS 100, 
piezosystem jena) focussed the excitation light and collected the fluorescence. The sample 
was scanned by a 2D piezo stage (PXY 80 D12, piezosystem jena). Collected fluorescence 
light was separated from the excitation laser by a dichroic mirror. Then, it was either 
spatially filtered and sent to a grating spectrometer or it was spectrally filtered by RG630 
glass, split by a 50/50 beam splitter, and directed towards two pinholes. These pinholes 
could be independently adjusted to spatially filter out different spots of the image plane, 
so that cross correlations could be measured. After the pinholes the light was directed to 
avalanche photodiodes (SPCM-AQR-14/SPCM-AQRH-14, Perkin Elmer). 

To perform correlation measurements at the waveguide ends, an oil-immersion microscope 
objective (UPlanSApo 60XO, Olympus) was used for collecting the light, while different air 
objectives were used to focus the excitation laser onto the sample from the other side. The 
air objectives were mounted on a 3D piezo system (PXY100 ID and MIPOS 100, piezosystem 
jena) to enable scanning of the excitation spot while the sample was kept at a fixed position. 
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Also, the 50/50 beam splitter within the HBT setup could be replaced by a D-shaped mirror 
in an image plane in cases where cross-correlation measurements with high count rates were 
needed. 

Q-factor measurements. The mode structure of the disc resonators was analysed by 
means of coupling in a tunable external-cavity diode laser (Velocity Series, New Focus) at 
around 770 nm via the evanescent fields of a tapered optical fibre. By tuning the frequency 
of the laser over the distinct whispering gallery modes, light was coupled into the resonator 
and different modes could be observed as Lorentzian shaped dips in the transmitted power. 
The polarization of the incoming light was chosen to maximize the coupling depth. For 
normalizing the data sets, we first performed a reference scan with an uncoupled fibre taper 
and compared the results to the situation when the same taper was coupled to the resonator. 
This provided direct access to the different coupling depths and cancelled out overlaying 
power modulations caused by slightly non-adiabatic taper transitions. 

The Q-factors were calculated from the dips by a Lorentzian fit function with an additional 
linear term to better match the local environment of the resonances. The fibre taper used 
had a waist diameter of 1.5 |xm. Measurements were performed in full contact. 
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SPECTRAL MODE ANALYSIS 



For a quantitative analysis of the spectrum measured from the background fluorescence 
by a grating spectrometer (see Figure [6]) we performed a segmented Fourier transformation 
on the data. We divided the set into 25 nm wide segments and performed numerical fast 
Fourier transformations on each of them individually. Afterwards we matched the frequency 
amplitude ranges by normalizing the different segments to each other, colour coded the 
resulting profiles, and merged the different wavelength ranges to obtain the plot shown in 
Figure [7j The different Fourier components of the resonator's whispering gallery modes 
(WGMs) are clearly visible as bands evenly spaced by an integer multiple of the lowest 
observable base frequency. This frequency component is directly related to 1/FSR and 
corresponds well with the theoretical estimated value based on the size of the discs. The 
black dotted lines are fits to this value and its different higher-order multiples. 
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FIG. 6. Spectrum of the background fluorescence when excited through one waveguide end and 
collected on the other (see main text) 
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FIG. 7. Normalised segmented Fourier transform showing the mode structure of the resonator's 
WGMs. The black dotted lines represent the different Fourier components of the collected spec- 
trum. They are fits to the expected Fourier frequencies of the resonator's modes and connected 
to the wavelength dependent FSR of the system. The lowest frequency component directly relates 
to 1/FSR. Within the determined spectral window a constant refractive index of around 1.5 was 
assumed for resonators with 20 \im diameter. 
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TUNING OF THE RESONATORS 



FIG. 8. Permanent tuning of a resonator under vacuum 



Figure [8] shows the tuning of a resonator under vacuum conditions. To tune the resonator, 
a 405 nm laser was focussed on its rim. The spectra shown were acquired from bottom to 
top in equal timesteps of 30 s. After some steps the mode begins to shift. It is worth noting 
that tuning of the resonator is possible even without oxygen present. 
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BACKGROUND CORRECTION 



To correct the measured g^(r) of an NV-centre for the collected background emitted by 
the resonator material, the following four steps were applied: 

1. Measurement of the normalised autocorrelation of the emitter with background 
a) and the normalised autocorrelation of the background gj?\r) 

(Fig. |b). 

2. Measurement of the background intensity relative to the emitter with background b. 
We measured b « |. 

3. The background corrected autocorrelation of the NV-centre is calculated according to 

4. The background corrected and normalised autocorrelation of the NV-centre is gf£j (r) = 
Iim^oo^(r)/^(t) (Fig. 0c). 

This background correction is exact, as calculated in the following. 

In the considered case of emission from a single NV-centre with background, the photon 
number operator can be decomposed into the sum of a part originating from the NV-centre 
n nv = a>nv°"nv an d the background contribution n& = a£a&. With the proper normalization 
g( 2 \oo) =< n >= 1 the of the emitter background system can easily be written down: 

g V)(r) = < [al(t) + a\{t)\[a{ v {t + r) + a\{t + r)\ 
x[a nv (t) + a b (t)][a nv (t + r) + a b (t + r)] > 
= < a) nv {t)a) nv {t + T)a nv {t)a nv {t + r) > 
+ < a\{t)a\(t + r)a b (t)a b (t + r) > 

+2 < ai v (t)a nv (t)al(t)a b (t) > (1) 

Now we make use of a simple mathematical trick called first order correlation expansion^ 
and define the correlation: 

< alajdiCLj > c := < a\a^aiaj > — < a\ai >< aja^ > . (2) 

For uncorrected operators ctj/j like the background operators the correlation vanishes and we 
get < alojdiaj >=< a\a,i >< OjCij >. Strictly, this is true only for the case of cw excitation. 
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£ +6 (r) (Fig. 



Nevertheless, as we will argue later on the procedure can also be used for measurements 
under pulsed excitation. Combining Es. 1 and Eq. 2 we find: 

^( 2 )( r ) = < a { v {t)a nv {t) > 2 + < a{ v {t)a{ v {t + r)a nv (t)a nv (t + r) > c 
+ < a\{t)a b (t) > 2 +2 < al v (t)a nv {t)a\(t)a b {t) > 
= (< a{ v {t)a nv {t) > + < a\{t)a b {t) >) 2 

+ < a{ v {t)a{ v {t + r)a nv (t)a nv (t + r) > c . (3) 

In the square root of the function the emission of the background n b is separated: 

\fgVHj) = < a) nv (t)a nv (t) > + < a\{t)a b {t) > + 

< al v (t)al v (t + r)a nv (t)a nv (t + r) > c 



2(< aL(t)a nv (t) > + < a\{t)a b {t) >) 
+ higher order terms {HOT) (4) 



= <n nv (t) > +b-yJg ( b 2 \T) + 

< al v (t)al v (t + T)a nv (t)a nv (t + r) > c 
' 2 
+ HOT, (5) 



t< 



where b is the background intensity relative to the total count rate and \J g^ (r) = 
al(t)ai,(t) > the autocorrelation of the Poissonian background. Here, it gets possible to 
correct the measured function for the independently measured g^ function of the back- 
ground. We can identify this background corrected function with the square root of the 
NV-centers normalised autocorrelation function: 

4WHr) = < n nv (t) > ^^H^ + ^aUVaUt + r) >< + R(yr 

At 

= 1 + < <4,(*) a L(* + T)a nv (t)a nv (t + r) > c HQT 

2 

= < a\iv{t)a\iv{t + r)a nv (t)a nv (t + r) > c 



= V < n nv (t) >< n nv (t + r) > + < aL(t)aL(t + r)a nv (t)a nv (t + r) > c 
= \J< n nv (t)n nv (t + t) > 



0") (6) 

The remaining function g( 2 ) can be is renormalized according to gf£} (r ) = lim^^ g^j (r ) / gf£} (t) . 

As mention, on the first view the above calculation seems to be valid only in case of mea- 
surements with cw excitation, since in pulsed excitation the background is not uncorrelated. 
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FIG. 9. Raw data used for the background correction, a and b are signal and background, 
respectively. Both were normalizes to an amplitude of one at long times, c shows the resulting 
antibunching. 

Anyhow, in pulsed excitation the measured can be decomposed into functions defined 
only at certain times after the repeated laser pulses: 

# (2) (r) if r = nAt + i;nGiV 
undefined else 

(2) 

For each g\ (r) the background is completely uncorrelated. The above derivation of the 
background correction is valid and gf$(r) can be reconstructed using the g^^r) functions. 



.(2) 
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